Ferro-Orbital Order and Strong Magnetic Anisotropy in the Parent Compounds of 

Iron-Pnictide Superconductors 
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The puzzling nature of magnetic and lattice phase transitions of iron pnictides is investigated via a 
first-principles Wannier function analysis of representative parent compound LaOFeAs. A rare ferro- 
orbital ordering is found to give rise to the recently observed highly anisotropic magnetic coupling, 
and drive the phase transitions — without resorting to widely employed frustration or nesting picture. 
The revealed necessity of the additional orbital physics leads to a correlated electronic structure 
fundamentally distinct from that of the cuprates. In particular, the strong coupling to the magnons 
advocates active roles of light orbitons in spin dynamics and electron pairing in iron pnictides. 

PACS numbers: 74.25.Jb, 74.25.Ha, 74.20.Mn, 71. 15. Mb 
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With a striking similarity to the cuprate families 
recently discovered high-temperature superconductivity 
in iron pnictides emerges upon suppression of magnetic 
order via additional charge doping [3, 0]. This new 
example of close proximity of superconductivity to the 
magnetic phase has raised again the long standing ques- 
tions concerning the intimate relationship between these 
two seemingly different (or even exclusive) phases of 
solids 0, i, E 0, i, S M, El, H E, 0, m • Following 
the cuprate research, a big thrust of current efforts has 
been to establish superconductivity in these compounds 
via magnetic correlations, despite the apparently diverse 
perspectives of the magnetism itself [3, S @| ■ Clearly, a 
solid understanding of the parent compounds and their 
magnetism is an essential first step toward a convincing 
resolution of superconductivity in the doped systems, es- 
pecially within the heavily discussed spin fluctuation sce- 
nario of pairing 0, i, S, 0, S, i, M, [Hi, [13, El, 0, El ■ 

The parent undopcd compounds of the iron pnictides 
have a quite unusual in-plane magnetic structure. Un- 
like the antiferromagnetic (AF) magnetic moments along 
both the X and y directions in the copper oxides (the G- 
AF structure), the magnetic moments in the iron pnic- 
tides are only AF in the x direction, but align ferromag- 
netically along the y direction. In the local picture, this 
curious stripe-like (or C-AF) structure is currently ex- 
plained [13, [3, [3] by requiring a strong next nearest 
neighbor (NNN) Heisenberg AF coupling, J2 , larger than 
half of the nearest neighbor (NN) AF coupling, Ji, such 
that the G-AF structure favored by Ji is suppressed en- 
ergetically. The competition between C-AF and G-AF 
implies a strong magnetic frustration 13, 13, 14 1. This 
frustration was argued to account for the observed small 
iron spin moment ~ 0.36 [3 and to promote super- 
conducting order via relief of the magnetic entropy [12| . 
In addition, a fluctuating electron nematic order was pre- 
dicted in the spin frustration case [l3| to account for the 
structural transition at slightly higher temperature [3], 



and was used to support the close relation between the 
physics of the cuprate and the iron-pnictide supercon- 
ductors [l3 |. 

However, if the system is really that frustrated, 
the rather high magnetic transition temperature 
(Tn K) [J^ would be difficult to understand. Fur- 

thermore, against the overall symmetry of the system, a 
surprisingly strong anisotropic NN coupling in the x and 
y direction was identified very recently from the inelastic 
neutron scattering measurements on CaFe2As2 \m, con- 
sistent with results of recent DFT calculations [3, [lit . 
This enormous anisotropy suggests that a strong rota- 
tional symmetry breaking has taken place prior to the 
magnetic ordering, setting a very stringent constraint 
to the correct microscopic understanding of the mag- 
netism. These results have been considered as direct 
evidence against [3] the above local Heisenberg picture 
and against [3] the alternative spin-density-wave pic- 
ture for itinerant electrons with "nesting" Fermi sur- 
face [3, [3, [1, [3, [l3l • Clearly, a comprehensive new picture 
is urgently needed to explain the microscopic origin of 
the magnetic structure together with the observed strong 
anisotropy, and in particular to address the potential ad- 
ditional symmetry breaking, considering its profound im- 
plications to the electronic structure and to the supercon- 
ductivity. 

In this Letter, we report a first-principles Wannier 



function analysis [16|, |17[ of the electronic structure of 



the representative parent compound, LaOFcAs. A purely 
electronic fcrro-orbital order is found to spontaneously 
break the rotational symmetry and drive the observed 
magnetic and structural transitions without resorting to 
the widely employed Fermi surface nesting or magnetic 
frustration. In great contrast to the cuprates, our study 
reveals the essential roles of the orbital degree of free- 
dom (especially its short-range correlation) in the elec- 
tronic structure of iron pnictides. In particular, the ex- 
pected light mass of orbiton and its strong coupling to 
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FIG. 1: (color online), a, Side view (left) and top view 
(right) for the Fel 3dxz (upper panels) and 3dyz (lower pan- 
els) Wannier orbitals obtained from nonmagnetic calculations, 
colored by their positive (red) and negative (blue) gradient, b, 
Sketches of top view of these orbitals on two nearest neighbor- 
ing Fe atoms along the x direction (Fel and Fe2). The solid 
(dashed) lines denote the tails above (below) the Fe plane. 



the magnon advocate its active roles in magnon decay 
and electron pairing. 

The electronic structure of LaOFeAs is calculated 
within local spin density approximation (LSDA) and its 
local interaction extension (LDA+U) of density func- 
tional theory, implemented via full potential, all-electron, 
linearized augmented plane wave basis [lit . The atomic 
positions are relaxed at the experimental lattice con- 
stants of the undistorted lattice at 175 K Q. A unit 
cell containing four Fe atoms in the plane is adopted to 
accommodate the observed C-AF structure. The x axis is 
chosen along the AF ordered direction of C-AF. A set of 
energy-resolved symmetry-respecting Wannier functions 
(WFs) is then constructed 3, 17 1 that spans the com- 
plete Hilbert space within 3eV of the Fermi level. 

In significant contrast to the maximally localized WF 
method [11] used in the repvious studie s \ldl . [20j , our con- 
struction exploits the "gauge freedom" [l9[ of the WFs to 
achieve localization within constraints of the point-group 



symmetry [16|, [ITf . Specifically, the cores of our Fe dxz 



and dyz WFs extend toward the directions of the mag- 
netic ordering, not the As atoms. Consequently, the one- 
particle density matrix, pij =<WF,;|p|WFj >, is locally 
diagonalized automatically, allowing direct detection of 
spontaneous orbital polarization 21 1 . 

Illustrated in Fig. [T^ are two of the resulting WFs most 
relevant to our further analysis (Fe 3dxz and 3dyz), as 
they are the only d-orbitals anisotropic in the xy direc- 
tions. To simplify the visualization, a sketch illustration 
of these WFs at the two neighboring sites is given in 
Fig. [T}d. The dramatic effects of hybridizing with tetra- 
hedral positioned As Ap orbitals can be clearly observed 
from the anti-bonding "tails" of the WFs. Notice, in 
particular, how such hybridization bends the tails of the 
WFs perpendicular to their original directions, and sig- 




FIG. 2: (color online). Electronic band structures of (a) non- 
magnetic and (b) C-AF configurations, represented in the 
nonmagnetic Brillouin zone. The weight of the Wannier Fe 
3da;2 and 3dyz orbitals are presented via blue and red colors, 
respectively. 



nificantly modifies the local point group symmetry. In 
great contrast to the normal two-dimensional character, 
the dxz and dyz WFs of LaOFeAs present a clasp-like 
shape around the iron center. Serious consequence of 
this hybridization will be discussed in detail below. 

With the help of these WFs, the LSDA band struc- 
tures for non-magnetic (NM) and C-AF LaOFeAs are 
compared in Fig. [21 To explicitly illustrate the effects 
of broken periodicity, the band structures are presented 
in the reciprocal space of the nominal NM unit cell con- 
taining only two Fe atoms. In this representation, addi- 
tional gap openings and "shadow bands" can be clearly 
observed in the C-AF case, whose intensity reflects the 
strength of these bands' coupling to the magnetic order. 

Notice that the bands near the Fermi level (zero en- 
ergy) primarily consist of the above mentioned Fe dxz 
and dyz WFs (c.f.: Fig. [1]). In the NM case, these two 
orbitals are degenerate as guaranteed by the point-group 
symmetry. In the C-AF case, on the other hand, a very 
large splitting in the broad vicinity of the Fermi level 
is observed involving the Fe dyz orbital, while the Fe 
dxz orbital is only weakly affected. This indicates that 
only one of the two orbitals is significantly involved in 
the magnetic coupling in the C-AF configuration, and the 
symmetry between the orbitals is broken. Indeed, as the 
magnetic order develops, the spin polarization (obtained 
from Pij) is found to be much stronger in the Fe dyz 
orbital ('^ 0.34/iB) than in the dxz orbital (~0.15 /zb), 
due to its loss of electron occupation in the spin minority 
channel. Introduction of a realistic moderate local inter- 
action {U=2 eV and Hund's coupling J= 0.5 eV) further 
enhances this effect and increases the Fe dyz polarization 
to 0.58 while the moment in the Fe dxz remains small 
(0.23 ^b) with both spin channels heavily populated. 

In addition, the gap opening, A, of the dyz orbital is 
found to be comparable to the band width, W, as shown 
in Fig. [2j Such a large gap is commonly encountered 
in strongly correlated systems where the magnetism is 
more conveniently described by interacting local mo- 
ments, rather than itinerant electrons with nested Fermi 
surface. Therefore, wc will proceed below to build a local 
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TABLE I: Onsite energy (first row) and hopping integrals 
among Fe 3d Wannier orbitals for the nonmagnetic case (in 
eV). Fe2 and Fe3 are the NN and NNN of Fel (c.f. Fig.[3|. 
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correlated picture for the electronic structure. 

Further insights into the microscopic local processes 
are revealed by transforming the DFT Hamiltonian of 
the NM case to the Wannier basis, as given in Table [D 
Unexpectedly, the leading hopping paths of Fe dyz orbital 
are to the neighboring dyz and d^-i.-yi orbitals along the x 
direction, perpendicular to its original direction. This is 
anti- intuitive since within the simple cubic symmetry the 
former would have been very weak "7r-bond"-like, while 
the latter would have been symmetry forbidden and iden- 
tical to zero. By contrast, the supposedly stronger "cr- 
bond"-like NN d^z — d^z hopping along the x direction is 
remarkably weak. These features are qualitatively differ- 
ent from those used in the previous studies [2^ . [isj . Our 
distinctly different results originate mathematically from 
the change of direction in the WFs' hybridization tails, 
as visualized above in Fig.[T] [^H. Physically, this reflects 
the dramatic influence of the tetrahedral positioning of 
As 4p orbitals on the Fe 3d orbitals, and reveals the im- 
portance of As atom positions and Fe-As phonon modes 
in the electronic structure in general. 

A simple "minimal" picture of the low-energy physics 
now emerges from the above analysis that elucidates 
the nature of C-AF magnetic structure together with 
the observed large anisotropy. Let's focus only on the 
Fe dxz and dyz orbitals, as they are the only d-orbitals 
anisotropic in the xy directions. As illustrated in Fig.[3l 
given almost doubly-occupied Fe dxz orbitals and almost 
singly-occupied, spin-polarized Fe dyz orbitals, the dyz 
orbitals prefer AF alignments along the directions of the 
efficient hopping, to benefit from the kinetic energy (the 
"super exchange" ) . Within the over-simplified strong 
coupling limit [2l| to the second order in the hopping 
parameters in TablelH the leading AF magnetic couplings 
among the dyz subspace are the NN coupling along the 
X direction, Ji^,, and the NNN couphng, J2 ~ 0.4 Ji^. 
In comparision, the NN coupling along the y direction, 
Jiy ~ 0.1 Jix, is insignificant. This large anisotropy 
is in good agreement with the current experimental [5| 
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FIG. 3: (color online), a Schematic of C-AF magnetic struc- 
ture with highly anisotropic NN coupling due to orbital or- 
dering, b, c The kinetic energy AE in the ferro-orbital and 
staggered orbital structures, t and U denote the hopping pa- 
rameter and intra-orbital Coulomb repulsion, while U' and 
Jh denote the inter-orbital repulsion and Hund's exchange, 
respectively. 



and theoretical observations [l5| . Clearly, the anisotropy 
owes its origin to the orbital degree of freedom, as the 
rotational symmetry breaking takes place via orbital po- 
larization. Therefore, modeling the magnetic structure 
with a standard Heisenberg model [1, 0, [l^, 14. 15| 
would suffer from its very limited applicability 2J] , as it 
lacks flexibility to break the rotational symmetry spon- 
taneously, or to adjust the magnetic coupling strength 
according to the orbital structure. 

As shown in Fig. [3l following these two leading AF 
couphngs Jix and J2, the observed C-AF structure is 
naturally established locally, without resorting to Fermi 
surface nesting. The seemingly ferromagnetic allignment 
along the y axis results primarily from the NN and NNN 
AF alignment across the columns. Contrary to previous 
theoretical explanation [l^ . [isl . [l^ . in our picture the 
observed C-AF structure is not frustrated or competing 
with any other magnetic structure {e.g.: G-AF), and thus 
can sustain a high transition temperature Q . 

In order to maximize the kinetic energy gain via the 
super- exchange processes, the orbitals not only have to 
be polarized, but also need to be ordered. (Other mag- 
netic/orbital configurations were found to have higher 
energies.) Indeed, in the above picture, all the sites are 
polarized the same way with the dyz orbital being less 
occupied and more spin polarized. This can be consid- 
ered an example of "ferro-orbital order" [i^, [2^. The 
formation of this rare orbital order can be understood by 
noting again the anti-intuitive NN hopping path along 
the X direction (c.f. Table |T|, dominated by only hop- 
ping between dyz orbitals without dxz-dyz cross hopping. 
As illustrated in Fig. ^p, the best way to utilize the ki- 
netic energy in this case is indeed the ferro-orbital, AF 
magnetic alignment, since one electron from both sites 
benefit from the kinetic energy. In comparison, the more 
common staggered-orbital, ferromagnetic alignment {e.g. 
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in undoped manganites) can utilize efficient hoppings in 
only one channel (c.f. Fig. [S};), despite the additional 



benefit from the intra-atomic interactions [2l|. That is, 
the unique hopping path leads to a rare ground state of 
the undoped iron pnictides consisting of coexisting, co- 
operative ferro-orbital and C-AF orders. 

The ferro-orbital phase is in fact directly responsible 
for the tetragonal-to-orthorhombic lattice transition ob- 
served at Ts ~155 K As soon as long-range ferro- 
orbital order takes place at Too{= Tg), dxz orbitals be- 
comes more occupied macroscopically (c.f. Fig. lead- 
ing to a longer bond in the x direction, in agreement with 
the experiment [26j . In general, it is very rare to have AF 
bond longer than the ferromagnetic one in the late tran- 
sition metal compounds, since the AF super exchange 
grows on almost half- filled orbitals (with less charge). In 
iron pnictides, this unusual long bond is realized only be- 
cause of the perpendicular extension of the hybridization 
tails discussed above. This subtle behavior distinguishes 
our result from a previous study [i^], in which the AF 
bond is less populated and thus would be shorter. 

It is important to further clarify the weak coupling 
strength between the elastic lattice distortion and the 
orbital order. When the lattice is set to frustrate the 
orbital order in our calculation by exchanging the lattice 
constants in the x and y directions, only less than 10 meV 
cost per Fe is found, an order of magnitude smaller than 
the super exchange energy. In essence, the orbital order 
originates almost entirely from the electronic energy, and 
the lattice simply follows the orbital order. On the other 
hand, such a weak coupling should allow a strong short- 
range orbital correlation even above Ts , offering a natural 
explanation to the large anomalous signal extending to 
40K above Ts in the thermal expansion measurement 27 1 . 

Since it is the same kinetic energy gain that drives 
both the magnetic and orbital orders, they arc thus 
strongly coupled to each other, as in the Kugcl-Khomskii 
model [25| . This is the natural reason for the close prox- 
imity of observed Tn and Ts- Unlike the manganites, 
where the orbital order is further stabilized by the large 
lattice-orbital coupling ~0.9 eV [17|, the weak couphng 
to the elastic lattice mode here does not help raising Ts 
noticeably from Tn. This suggests a light effective mass 
of the orbiton, comparable to the magnon mass. Thus, 
an intimate relationship between the excitations in the 
spin and orbital channel is expected. Indeed, an efficient 
orbiton-assistcd decay offers a very interesting possible 
explanation to the observed large line width of magnon 
at large momentum Furthermore, the controlling role 
of Fe-As-Fe positioning in the unconventional hopping 
implies large isotope effects in both magnetic [28[ and 
orbital transitions. 

Our study also has several implications on the high- 
temperature superconductivity of doped iron pnictides. 
In great contrast to the cuprates, our results reveal the 
essential roles of the orbital degree of freedom in iron 



pnictides. Upon doping, the long-range magnetic and or- 
bital orders would naturally perish rapidly through the 
disruption of the above superexchange process by doping- 
induced charge fluctuation. (Recall that the above pro- 
cess maximizes with three electrons in the xz / yz complex 
in average.) Nonetheless, the short-range orbital correla- 
tion should persist deep into the underdoped supercon- 
ducting regime. In addition, the light mass of orbiton and 
its strong coupling to the magnon make orbiton another 
interesting participant in the pairing mechanism. Fur- 
thermore, the correlated nature suggested by our study 
indicates a much stronger local electron-boson coupling 
than the current mean-field estimation [2^ . Finally, our 
correlated picture also supports a strong coupling na- 
ture of superconductivity (with strong phase fluctuation) 
in the underdoped regime, in agreement with the recent 
observation of low super-fluid density [sot with relatively 
high T,. 

In summary, with our flrst-principles Wannier func- 
tion analysis, a rare ferro-orbital order is identified in 
LaOFeAs that breaks the rotational symmetry and gives 
rise to the recently observed giant anisotropy of mag- 
netic coupling in undoped iron pnictides. The observed 
magnetic and lattice structures result locally — without 
widely applied frustration or the need for Fermi surface 
nesting. In great contrast to the cuprates, our study re- 
veals the essential roles of the orbital degree of freedom in 
iron pnictides in general, and suggests active roles of light 
orbiton in magnon decay and electron-pairing. Finally, 
the pure electronic origin and the obtained large cou- 
pling strength advocate a correlated nature of undoped 
iron pnictides, and support the notion of strong coupling 
superconductivity in the underdoped regime. 

This work was supported by the U.S. Department of 
Energy, Office of Basic Energy Science, under Contract 
No. DE-AC02-98CH10886, and DOE-CMSN. 
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